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A SYSTEM AND METHOD FOR ENHANCING 
THE ACCURACY OF A LOCATION ESTIMATE 

CROSS REFERENCES 

[0001] The present application is co-pending with and claims priority benefit of 
provisional ^plication entitled "Geolocation of Mobile Appliances", Appl. S.N. 
60/418,342 and filed on October 16, 2002, the entirety of which is hereby incorporated 
h^in by reference. 

[0002] The present application is related to and concurrently filed with 
applications titled "A NETWORFC OVERLAY GEO-LOCATION SYSTEM WITH 
SMART ANTENNAS AND METHOD OF OPERATION" SN #, "WIRELESS 
COMMUNICATION NETWORK MEASUREMENT DATA COLLECTION USING 
INFRASTRUCTURE OVERLAY-BASED HANDSET LOCATION SYSTEMS" SN #, 
"NETWORK OVERLAY LOCATION SYSTEM AND METHOD FOR AIR 
INTERFACE WITH FREQUENCY HOPPING" SN #, "A SYSTEM AND METHOD 
FOR ESTIMATING THE MULTI-PATH DELAYS IN A SIGNAL USING A 
SPATIALLY BLIND ANTENNA ARRAY, SN#, and "SYSTEM AND METHOD FOR 
OPERATING A NETWORK OVERLAY GEO-LOCATION SYSTEM WTIH 
REPEATERS" SN #, filed October 16, 2003, the entirety of each of these applications is 
incorporated h^in by reference. 

BACKGROUND 

[0003] In a typical location estimation scenario, several receiving stations 
simultaneously derive measur^ients on the emitter signal, the emitter being a wireless 
transmitter, a mobile appliance such as a mobile phone. Personal Digital Assistant 
("PDA"), orpersonal computer with wireless cj5)ablitity. Ideally, any one of these 
measurements is a fimction only of the emitter location P , the receiving station locations 
Sj (where the subscripty denotes die station), and flie antenna configurations at the 
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receiving stations. A given receiving station may attempt to derive the bearing or Angle 
of Arrival ("AOA") at which the emitter is located. A different receiving station may 
attCTipt to derive the Time Of Arrival CTOA") of the emitter's signal. Some receiving 
station pairs may attempt to derive the Time Difference Of Arrival ("TDOA") of the 
emitter signal between the station pair. Other receiving station pairs may attempt to 
compute the Frequency Difference Of Arrival ("FDOA") of the emitter signal between 
the station pair. The form of the measurements is not restricted to the above; certain 
receiving stations may derive a multiplicity of the measurements indicated or other more 
exotic measurements. 

[0004] Most prior art location estimating systems have more receiving stations in 
place than minimally required. For example, if all receiving stations were to use TOAs 
to determine the emitter position, three stations would suffice in the ideal case. For a 
perfect estimate in this scenario, the three TOAs will be exactly correct (perfect) 
measurements. For receiving station using AOA, two stations would suffice in the ideal 
case. Where the number of receiving stations used in the estimate is above the ideal 
number, the estimate becomes an overdetemiined solution. The fact that such perfect 
measurements are never available in practice necessitates the use of an excess of 
receiving stations and a resulting excess of measurements and an overdetermined 
solution. Given the excess measurements, a location estimate is derived to best fit the 
measurements. 

[0005] In a mixed mode system, a combination of TOAs, TDOAs, AO As, 
FDOAs and other measurements are combined to estimate the emitter position. Whether 
mixed mode or not, the same principle of using (or attempting to use) an excess of 
receiving stations is appUed to generate a more reliable estimate of the emitter location. 

[0006] In most prior art systems, the mathematical approach taken to derive a 
location from such an excess of measurements assumes that each of the measurements is 
the perfect measurement corrapted with Guassian noise widi some known statistics. A 
more detailed explanation of this approach can be found in M.Wax, "Position location 
from sensors with position xmcertainty", IEEE Trans. Aero., Elect. Syst. AES-IP, no. 2 
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(Sept. 1983), 658-662; D.L Torrieri. "Statistical Theory of Passive Location Systems", 
IEEE Trans. Aerosp. Electron. Syst. AES-20, no. 2 (Mar. 1984), 183-198; Y.T. Chan and 
K.C. Ho, "A simple and efficient estimator for hyperbolic location", IEEE Trans. Signal 
Proc. 42, no, 8 (Aug. 1994), 1905-1915; W.H. Foy. "Position location solutions by 
Taylor series estimation", IEEE trans Aerosp. Electron. System AES-12, no. 2 (Mar. 
1976), 187-194; R.a Stansfield, "Statistical theory of DP fixing", Joum. lEE 94, part 
Ilia (Oct. 1947), 762-770; the entirety of each is herein incorporated by reference. This 
technique has a long and established history and serves as the bedrock of location 
estimation. 

[0007] In actual systems, modeling the receiver station measurement as a perfect 
signal corrapted by noise is accurate only in a small minority of cases. The reason for 
this is tihat measurements (of any of the types indicated earlier) typically have biases 
which are rarely (if ever) reflected in the noise statistics. This measurement bias may 
stemfiomavariety of factors. One source of bias is instrumentation error. Another 
source of bias is signal multipath where a delayed signal, a reflection of the original 
signal, masquerades as the desired signal. 

[0008] In determining a location estimate, it is convenient to associate a weight 
with each measurement. This translates mathematically to either individual weights or a 
matrix that expresses inter-relationships among the measurements. Given that the biases 

are unknown and may (in the case of multipath) be functions of the knitter location P , a 
direct ma&ematical solution using both the initial weights and unknown biases is 
impossible: there are an infinite number of solutions. On the other hand, all known 
solution techniques that ignore tiie measurement bias result in an estimate that is itself 
biased. Since the location estimate is derived by a mathematical weighting of each 
measurement, the weights applied to each measurement have a strong effect on the error 
in the location estimate. Ideally, the weights applied to good measurements should 
always be larger than those applied to biased measurements. 

[0009] Several prior art attempts have been made to address the issue of biased 
measurements and are described in detail in M.P. Wylie and J. Houtzman, "The non-line 
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of sight problem in mobile location estimation". Proc IEEE S^Iintemattonal Conf. on 
Universal Personal Communications, vol 2 (Oct 1996), 827-831; L.Cong and W.Xuang, 
''Non-Line-of-Sight Error Mitigation in TDOA mobile location" Proc. IEEE Global 
Telecommunications conference voll (2001), 680-684; P.C. Chen, "A non-Iine-of-sight 
error mitigation algorithm in location estimating" Proc. IEEE Conf, on wireless 
Communications Networking, vol 1 (1999), 316-320; andNJ. Thomas, D.G.M. 
Cruickshank and D.LLaiirenson, "Performance of a TDOA-AOA hybrid mobile location 
system" 3G Mobile Communication Technologies Conf. Proc. 1 (Msa:, 2001), 216-220, 
all of which are incorporated herein by reference. These references describe approaches 
that identify the offending measurements and then either eliminate such measurements or 
model the offending measurements with a distribution different from the traditional 
approach. Some of these methods additionally require a large number of samples of a 
particular measurement to create a time-history of the measurem^t. 

[0010] A major problem with the prior art techniques is &at in practical systems 
die differentiation between biased and non-biased measurements is never clear-cut. 
Unlike in purely academic simulations, real life data reveals a continuum ranging from 
near perfect measurements to measuremaits with large biases. E:q)erimentation with 
schemes that attempt to isolate a particular biased nieasurement have shown that such 
schemes rarely work. These approaches have great difficulty in identifying the offending 
measurements when more than one receiving station is in error. 

[0011] For purposes of this disclosure, domiaant measurements are those 
measurements that most strongly influence the location estimate. When the dominant 
measurements have a smaller bias than the remaining measurements, the estimate 
generated by any of the traditional prior art solution techniques mraitioned previously is a 
better estimate than the estimate generated by the measurements with larger bias. That 
is, the estimate is far from perfect, but is skewed to favor the dominant measurements 
that within this scenario are desoibed as having a small bias. In practice, and especially 
when a fair excess of measurements is available, the dominant measurements actually 
have a smaller bias. 
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[0012] One reason for this is that the set of low bias measurements have a greater 
degree of self-coherence (less variation or greater mutual agreement) with respect to the 
mattiematics that generates the estimate; hence, tiiey are more likely to dominate the 
estimate. The skew in the estimate makes the non-dominant measurements have a larger 
offset with respect to the estimate than the dominant measurements. It is important to 
note that one or more non-dominant measurements may have a large weight The self- 
coherence of several lower weight measurements can thus dominate the location 
estimate, over-riding the non-dominant high weight measurement that is possibly in 
error. 

[0013] The disclosed subject matter capitalizes on this fact to adjust the weights 
applied to the measurements by the offsets from the hypothetical measurements, with the 
assumption that at each iteration the location estimate is exact. Thus, the relative weights 
applied to the dominant measurements are increased with respect to the non-dominant 
measurements. This process is recursively refined to generate further improvements in 
the location estimate. The dominant measurements may or may not be a majority of the 
measurements since only a few high weight measurements may dominate tiie estim^e. 
Conversely several lower weight measurements may prove dominant in terms of the 
location estimate, and hence their weights may increase in the next iteration of the 
process, while the weigjits applied to the other measurements may decrease. 

[0014] The method disclosed makes no attempt to determine the offending 
measurements. The goal, rather, is to improve the location estimate. However, the 
method can be used to identify tiie biased measurements subsequent to refining the 
location estimate. 

[0015] Therefore, it is an object of the present subject matter to obviate the 
deficiencies of the prior art and present a novel method and system for recursively 
refining location estimates by accounting for receiver bias where an overdetermined 
solution exists. 

[0016] It is also an object of the present subject matter to present an improved 
method for refining a geo-Iocation estimate of a wireless transmitter emitting a signal that 
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is received by a predetennined number of sensors that is greater than the minimum 
number of sensors required to obtain the geo-location estimate. The improvraient may 
include incorporating bias error in the signals received at the sensors and updating the 
geo-location estimate by recursive analysis of the bias error to thereby refine the geo- 
location estimate. 

[0017] It is still an object of the present subject matter to present a method for 
estimating the geo-location of a wireless transmitter emitting a signal that is received by 
a plurality of sensors in a geo-location system which includes a geo-location estimation 
device which provides an overdetermined geo-location solution for the wireless 
transmitter. The method includes measuring an attribute of the emitted signal to thereby 
create a sensor signal at the sensor and sending the sensor signal to the geo-location 
estimation device. The method may also include receiving the plural sensor signals, 
associating with each sensor signal a separate initial predetermined weight value to 
thereby provide a plurality of initial estimation signals, determining an initial estimate of 
the geo-location of the wireless transmitter from the initial estimation signals, and 
modifying the weight value associated with the dominant sensor signals relative to the 
weight value associated with the non-dominant sensor signals to thereby provide a 
plurality of refined estimation signals. 

[0018] It is another object of the present subject matter to present a novel metiiod 
for estimating the geo-location of a wireless transmitt^ emittmg a signal that is received 
by a plurality of sensors in a geo-location system which includes a geo-location 
estimation device which provides an overdetermined geo-location solution for the 
wireless transmitter as a Amotion of sensor signals determined firom an attribute of the 
received signal at the plurality of sensors. 

[0019] It is yet another object of the present subject matter to present a novel 
system for estimating the geo-location of a wireless transmitter emitting a signal that is 
received by a plurality of sensors in a geo-location system which fiirther includes a geo- 
location estimation device which provides an overdetermined geo-location solution for 
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fhe wireless transmitter as a function of sensor signals determined firom an attribute of the 
received signal at the plurality of sensors. 

[0020] These and other advantages of the disclosed subject matter over the prior 
art will be readily apparent to one skilled in the art to which the disclosure pertains &om 
a perusal of the claims, the upended drawings, and the following detailed description of 
the preferred embodiments. 

BRIEF DESCRIPTON OF THE DRAWINGS 

[0021] Figure 1 is an illustration of the elements used to estimate a location in a 
TO A geolocation system according to an embodiment of the disclosed subject matter. 

[0022] Figure 2 is an illustration of the elements used to estimate a location in a 
AOA geolocation system according to an embodiment of the disclosed subject matter. 

[0023] Figure 3 is a flow chart representing an embodiment of the disclosed 
subject matter. 

DETAILED DESCRIPTION 

[0024] A mathemetical description follows to aid in the description of the current 
subject matter. The set of measxjrements from the receiving stations is denoted by m . 
These measurements are determined from an attribute of the received signal from the 
mobile emitter. The receiving stations are typically sensors co-located with base stations 
in a wireless conmiunication system. The sensors can also be located at repeater stations 
within the communication system or remotely located. The elements of this colimm 
vector m are the individual measurements, possibly of several different types, such as 
TOA, TDOA or AOA, FDOA etc. Each component in the vector is associated with one 
or more of the receiving stations 7, where J denotes (indexes) the station. The self and 
inter-relationships between the measurements are denoted by the weight matrix Wjt 
where k denotes the iteration number, k is equal to unity at the outset and increases in 
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steps of 1 at each iteration. The matrix formulation is adopted for simplicity since the 
actual inter-relationships between the measurements could have a more complex form* 
In the simplest case, the weight matrix W^t is a diagonal matrix expressing the variance 
associated with individual measurements. 

[0025] The emitter location estimate is denoted by , and the receivmg station 
locations are denoted by the column vector s . The emitter location estimate is given as: 

[0026] Pk-^fis, m. Wit) ;^ =1,2,3,.... (1) 

Thus the location estimate is a function of receiving station locations, measurements and 
weights. The subscript k denotes the iteration number. 

[0027] Here, / (? , iw , Wjt.) is the function used to generate a location estimate 
given the set of measurements and their self and inter-relationships. Depending on the 
approach talcen, f{s, m , Wjk) could be one of many diffCTcnt functions. For example, 
any one of the functions expressing the location estimate in terms of the sensor positions 
and the measurements as described in the prior art may be used. 

[0028] It is clear tiiat the established methods encapsulated by equation (1) need 
to use a particular weight matrix. In our formulation, the initial weight matrix Wi, is a 
known weight matrix which can be predetermined based on prior Imowledge or attributes 
of the received signals such as the Signal to Noise Ratio (SNR), or even in the most 
trivial case, an identity matrix. The initial location (P^) may be derived from an expUcit 
solution such as in the prior art solutions discussed previously or derived using a specific 
algorithm satisfactory to the practitioner of this art. 

[0029] Methods embodying the disclosed subject matter iteratively update or 
modify the weight matrix in accordance with: 

[0030] Wjtfi =g(W/t, 5, m,P4);lc =1,2,3,... (2) 
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where flie function gfWjfc, s^nty ) niodifies the weight naatrix depending on the 
agreement between the individual measurements in m , the receiving station locations Sj 

and the emitter location estimate at the ^ step, . 

[0031] Equation (2) summarizes the general approach embodied by the disclosed 
subject matter. To provide more detail for a particular case, an embodiment of a location 
systCTi using only TOAs is illustrated. 

[0032] Figure 1 is a representation of a geolocation system using TOA, showing 
the respective elements of equation (1). The receiving stations are designated as 101- 

104, the actual mobile emitter location P is shown as 1 10 while the estimated mobile 

emitter location at recursion or iteration ^ , is shown as 1 1 1 . 

[0033] Let the sensor locations be denoted by Sy, where y= 1, 2...N; iV denoting 
the niunber of sensors; in Figure 1, N=4. 

[0034] An embodiment of equation (2) for this case is then: 

[0035] Witf/ = Wife + a ( - - (3) 

[0036] in which 

[0037] <ic K^l,k h'k — ^V/ifc)^ = (^l»k fe-k — Uk)^ 

[0038] where *T*' denotes matrix transposition, 

and a is a constant matrix which in its simplest form is a constant scalar times 
an identity matrix and can be theoretically or empirically determined or provided by the 
geo-location system operator as a fixed setting or in real time. 

[0039] As a second illustration of the method, consider a scheme that uses only 
TDOAs. In this case a possible embodiment of equation (2) is 
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[0040] r^^, =W, +aim-r,X^-Tj 



(4) 



[0041] where each element in m denotes the TDOA measurements with respect 
to a particular receiving station pair. 



[0042] ^-jfc'Ht* 



[0043] the (AT - i) X (N) matrix H is given by: 



[0044] 



10 0 
0 1 0 



0 -1 
0 -1 



0 0 0 



1 -1 



and the same reference station (station JV^in fbia case) is used in the formation of &e 
TDOA measurements. 

[0045] A third embodiment of equation (2) is shown in Figure 2, where the 
location estimation scheme uses AOA. In this case a possible embodiment of equation 
(2) is 



[0046] jn*. = +<xi^-nX^-rJ 



(5) 



[0047] where each element in m denotes the AOA measurement with respect to 
a particular receiving station, and 



=tan 



-1 



where subsoipts x and y indicate the individual components along the x-axis and y-axis in 
a 2-dimensional coordinate scheme y j indexes the station and k indexes the iteration. 
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[0048] It will be understood by those of skill in the art that although the above- 
described embodiments describe systems with four sensors, systems with other than four 
sensors are contemplated by the present subject matter consistent with the decription 
herein. 

[0049] Figure 3 shows a flow chart illustrating a typical implementation of the 
disclosed subject matter. Measurements or sensor data, where the latter is the precursor 
of the former along with an initial Weight matrix are supplied to block 3 10 where any of 
the prior art methods consistent with equation (1) are used tq produce an intial fix. The 

intial fix or location estimate Pj along with the measurements or sensor data are supplied 
to block 320 where measurement of&ets are determined. These offsets are used to 
update the weight matrix according to equation (2) in block 330. A new location 

estimate is generated with the updated or modified Weight matrix and the 

measurement m in block 340. The iterations continue until either Ihe Position estimates 
converge as shown in block 350, or altematively the iterations cease when some attribute 
of the Weight matrix converges. 

[0050] The iteration process can be implemented at the geo-location estimation 
device or in a subsystem in operational conununication with the geo-location device. 
The geo-location device is generally a processor that performs tiie mathematical 
manipulation of the signals fi-om the sensors. The mathematical manipulation can be 
accomplished with hardware and /or software. 

[0051] The technique of the current subject matter iteratively modifies the weight 
matrix W* as CT:pressed by equations (1) and (2) is functional over varying terrain, 
different receiving station configurations, and over widely varying cellular protocols. 

[0052] While preferred embodiments of the present inventive system and method 
have been described, it is to be imderstood that the embodiments described are 
illustrative only and that the scope of the embodiments of the present inventive system 
and method is to be defined solely by the appended claims when accorded a full range of 
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equivalence, maay variations and modifications naturally occurring to those of skill in 
the art from a perusal hereof. 
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